Test of molecular mode coupling theory: A first resume 
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We report recent progress on the test of mode coupling theory for molecular liquids (MMCT) for 
molecules of arbitrary shape. The MMCT equations in the long time limit are solved for supercooled 
water including all molecular degrees of freedom. In contrast to our earlier treatment of water 
as a linear molecule, we find that the glass transition temperature T c is overestimated by the 
theory as was found in the case of simple liquids. The nonergodicity parameters are calculated 
from the "full" set of MMCT-equations truncated at l co — 2. These results are compared (i) 
with the nonergodicity parameters from MMCT with l co = 2 in the "dipole" approximation n = 
n' — and the diagonalization approximation n — n' = 0, I — I' and (ii) with the corresponding 
results from a MD-simulation. This work supports the possibility that a reduction to the most 
prominent correlators may constitute a valid approximation for solving the MMCT equations for 
rigid molecules. 

PACS numbers: 61.25.Em, 64.70.Pf, 61.43.Fs, 61.20.Ja 
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I. INTRODUCTION 

The mode coupling theory (MCT) for supercooled sim- 
ple liquids proposed by Bcngtzelius, Gotze and Sjolander 
|jj interprets the glass transition as a dynamical transi- 
tion. This picture has been supported by many exper- 
iments on several glass formers (see e.g. || and Refs. 
therein) and more recently by detailed analysis of com- 
puter simulations for Lennard- Jones systems (see e.g. || 
and references therein). The signature of the dynamical 
transition, i.e. the asymptotical power laws, have been 
discovered also via neutron scattering, light scattering, 
dielectric relaxation and NMR in molecular glassformers 
(see H and references therein), stimulating the extension 
of MCT to molecular liquids. 

Two approaches for such an extension have been pro- 
posed recently for rigid molecules. Chong and Hirata 
j| introduced a theory based on a site-site description 
of the molecules. Their approach offers the advantage 
to be closely related to neutron scattering experiments. 
The structural information — which is a necessary in- 
put of the theory — can be readily obtained from theo- 
ries of molecular liquids able to predict partial structure 
factors, like the RISM || approximation. The second 
approach is based on the expansion of the orientational 
density into a complete set of functions, in analogy to 
the Fourier expansion of the density related to the trans- 
lational degrees of freedom. To distinguish the second 
approach from the site-site theory, it is called molecular 
mode coupling theory (MMCT). MMCT was derived for 
a single linear molecule in a simple liquid for liquids 
of linear molecules (for some application see |s) ) and 
for molecules of arbitrary shape [p|flO|] . MMCT allows to 
calculate the glassy dynamics for all orientational degrees 



of freedom, and it is closely connected to dielectric and 
NMR experiments. As reorientational motion is also very 
important in light scattering experiments 1 1 1 1^] it may 
also be helpful in their interpretation. The connection 
of MMCT to neutron scattering experiments has been 
discussed recently Q . 

The fundamental MMCT quantities are the time- 
dependent correlation functions 

Si n ,i' n '(q,m,t) = (p* ln (q,m,t)p Vn ,{q,m)) (1) 

of the tensorial density modes 

N 

Pln (q,m,t) = i\2l + 1)4 ^e^f'^^fi)). (2) 
i=i 

Here it is q = (0, 0, q) and I runs over all positive integers 
including zero, m and n take integer values between —I 
and I, V denotes the Wigner functions Jl3]. The reader 
should note that the correlators ([l]) are diagonal in m 
for q = (0,0, q). The MMCT equations of motion for 
the Laplace transform S(q, m, z ) = i S(q, m, t)e lzt , 
(Im z > 0) have been presented in a preceding paper. 
Here we focus on the unnormalized molecular nonergod- 
icity parameters 



F(q,m) = lim S(q,m,i) 

t — >oo 



lim zS(q, m, z) (3) 



which we calculated using the following set of equations 

F(q,m)= (4) 
[S-% m) + S- 1 ^, m)K(q, m)S" 1 (g, m)] _1 



1 



Kln,l'n'(<l,m) = (5) 
C(9)^'*(9) (sT 1 (<?)) ^n" w (6) 



" l lmn,l'm'n'\'J-! ~ (87T 2 ) 



oo <2+<?i 
|8-«i| 



H v rrC; 1 n 1 ,i 2 ™ 2 (wi92;mmim 2 ) 

rn 1 m 2 li'^ H»J 

^ 2 n 2) i 2 ^(92,m 2 ) (7) 

The memory function matrix m in Eq.(Q) represent the 
mode coupling approximation for the correlation func- 
tion of fluctuating forces. The index a labels the transla- 
tional (a = T) and rotational (a = R) currents, each of 
them consisting of three (spherical) vector components 
p e { — 1,0,1}. The vertex functions v are determined 
only by the matrix of the static molecular structure fac- 
tors S(q, m) and the number density p . Their explicit 
form has been given in Ref. 0. The coefficients q^(q) 
appearing in Eq.(||) are 



a = T, p = ± 1 
a = T, p = 
a = R, p ± 1 
a = R, p = 



(8) 



Since qj^ 1 — (due to the choice of q = (0,0,(7)) the 
transversal translational components (a = T , p = ±1) 
of the memory functions enter only indirectly (via the 
inversion of m(g)) and thus shall be neglected in the fol- 
lowing. 

The given set of equations includes all in- 

teractions between translational and rotational degrees 
of freedom in molecular liquids and thus this set is 
rather involved. Obviously its numerical solution poses 
a formidable task. 

As model system for our analysis we have chosen 
SPC/E water |Q. The water molecule possesses a 
twofold rotational symmetry (C2 V ) around the axis given 
by its dipole moment which has been chosen as the z-axis 
of the body fixed frame. 





MD 


MMCT diag 


MMCT dipole 


MMCT "full" 


T c 


200 K 


206 K 


208 K 


279 K 



TABLE I. Comparison of the critical temperatures T c for 
the MD-simulation and the theoretical calculations in diag- 
onal-dipole approximation, dipole-approximation and for the 
"full" theory. 



The z-axis and the x-axis define the plane which is 
spanned by the molecule. As discussed in detail in Ref. 
]j~0t the C2t,-symmetry can be used to simplify the equa- 
tions of motion by the restriction that n and n' are even. 
Prccccding publications on long time MD-simulations 
p6[ for this strong glass former showed that the centre of 
mass dynamics is in good agreement with the predictions 
of the asymptotic laws of MCT. Recent work [jl7| demon- 
strated that the signature of the dynamic transition can 
also be observed for the orientational degrees of freedom 
of the molecule. 

A first approach to solve Eqs. (|5|)-(0) has been given 
in Ref. Jl(|. Apart from the necessary truncation of 
the range of I by a cutoff l co for which we have chosen 
l co = 2 we introduced in fic[| two more approximations. 
In the dipole approximation, water was treated as a linear 
molecule oriented in the direction of its dipole moment. 
This approximation corresponds to set in Eqs. @-(0) the 
angular indices n, n' and the corresponding summation 
indices n i ,n' i to zero. In Ref. jl0| we have also studied 
the stronger diagonal- dipole approximation, which makes 
the further assumption that structure factors, nonergod- 
icity parameters and memory functions are diagonal with 
respect to the angular indices / and It is the main 
purpose of the present communication to present results 
obtained by solving the "full" set of MMCT-equations 
up to l co = 2 without any further approximation. The 
main difference to the approach in ref. jl(^] is that we take 
the nonlinear character of the water molecule serious. In 
particular, this means that the rotational motion of both 
protons around the dipolar axis is taken into account. 
Therefore, we can study the influence of this degree of 
freedom on the ideal glass transition temperature T c and 
on the nonergodicity parameters i 7 io,ro('?i fri) obtained in 
ref. |10|. In addition we also obtain the new parameters 
Fin t i' n i{q,m) with n and/or n' different from zero. 



II. RESULTS 

We have solved Eqs. (|[)-(0) iteratively on a grid of 100 
equispaced wave-vectors, up to 110.7 nmT 1 . One com- 
plete iteration requires 4 days of cpu time on a 533 MHz 
alpha workstation. The entire calculation to locate the 
critical temperature and the corresponding nonergodic- 
ity parameter, with a tolerance of 2.5 x 10~ 4 per point, 
requested more than 250 iterations. On a dedicated 4- 
nodes parallel machine it required about 250 days. The 
number of iteration at the critical temperature was 54. 
We found that at T = 2S2K the MMCT equations pre- 
dict a liquid phase, while at T — 272K unambiguously 
a glassy state is predicted. Within the chosen tolerance, 
we locate T C MCT at T C MCT = 279K. 

Table [| summarizes the critical temperatures at which 
a transition from ergodic to nonergodic behavior is found. 
While T c is almost equal for both approximations used in 



2 



ref. |Tc}] and quite close to the result of the MD-simulation 
we find that the critical temperature is overestimated by 
almost 50% using the "full" theory. Thus we confirm 
our supposition in ]l0| that the agreement for T c be- 
tween simulation and both approximation schemes was 
fortuitous. Such finding could also be expected from 
the static structure factors used as input. Besides the 
dominating diagonal structure factors Soo,oo(9,to = 0) 
and Sio.io (<?, to) the most prominent peaks are displayed 
by £io,2±2(<7, to) and S 2 ±2,2±2(q, to) (See Figures in Ref. 
p8|), which were neglected in [ ]l0| |. The overestimation 
of the critical temperature is common to the MCT for 
simple liquids J3| and seems to be a general deficiency 
of the mode coupling approximation. Thus we see that 
- although the overestimation of the critical tempera- 
ture is not desired — it is necessary to include all static 
structure factors with large amplitude to get a concise 
description in the MMCT framework. 
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FIG. 1. Critical nonergodicity parameters -£00,00(9, m = 0) 
for the centre of mass. The results from the MD simulation 
(symbols) are compared with the theoretical predictions of 
MMCT in diagonal-dipole approximation (short dashed line) , 
in dipole approximation (solid line) and the results obtained 
using the "full" set of MMCT-equations (long dashed line). 

Figure [I] shows the comparison of the unnormalized 
critical nonergodicity parameters -£00,00(9, to = 0) for the 
centre of mass correlations. The oscillations of the MD- 
result are captured well by all of our MMCT-results. As 
for the critical temperature we can also observe for the 
critical nonergodicity parameters that the dipole approx- 
imation and the diagonal-dipole approximation lead to 
nearly the same result. In the vicinity of the maximum 
of the structure factor the agreement between theory and 
simulation is improved by removing the additional ap- 
proximations but for the prepeak as well as for large 
q both approximation schemes perform better than the 
" full" theory. In the region of the prepeak and especially 



for the minimum between prepeak and main peak we ob- 
serve that the oscillations are less pronounced and the 
peak positions are slightly shifted. Exactly the same be- 
haviour can be found by a comparison of the static struc- 
ture factors at the different critical temperatures. Thus 
the worse performance of the "full" theory in compari- 
son with the approximation schemes can be attributed 
at least partly to the overestimation of the critical tem- 
perature because of which the static input of the calcu- 
lations does not reflect accurately the static structure at 
the "true" T c of the simulation. Apart from that, it also 
has to be taken into account that in spite of the com- 
putational effort we have made, the fixed point of Eq. 
(H)"(@) can n °t De determined with very high precision. 
Therefore, the result for the critical nonergodicity pa- 
rameters after 54 iterations may overestimate the exact 
one by a few percent. The latter reason may also be 
responsible for the worse performance of the " full" the- 
ory at the minima of Fbo,oo(<?, to = 0) as Nauroth |lj| 
found that the convergence of the iteration is extremely 
slow in these parts. Finally we want to mention that the 
disagreement at large q between simulation and all the- 
oretical calculations was also found for simple liquids || 
and was considered as a shortcoming of the mode cou- 
pling approximation. 
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FIG. 2. Diagonal critical nonergodicity parameters 
Fio,io(<Z, m) as calculated from the MD simulation (symbols) 
compared with the theoretical predictions of MMCT in diag- 
onal-dipole approximation (short dashed line), in dipole ap- 
proximation (solid line) and the results obtained using the 
"full" set of MMCT-equations (long dashed line). 

As in the case of the centre of mass correlators the q- 
dependence of the nonergodicity parameters Fio,io(9, to) 
(see fig. ^ is reproduced well by the theory. Here even in 
the vicinity of the maximum the approximations perform 
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better than the "full" theory. 

Figure ^ shows the comparison for the critical noner- 
godicity parameters -F 2 o,2o(<2S m )- Apart from the region 
of large wave vectors the " full" theory shows better agree- 
ment with the simulation than the two approximations. 
Thus one observes that the approximations, which ne- 
glect terms with n ^ and consequently I > 2, have 
stronger effect on the I = 2-correlators than on those for 
I = 1 or I = 0. Further one observes that the agree- 
ment between theory and simulation is less good for the 
I = 2-correlators than for those with I = or I = 1. The 
reader should note that the good agreement at large q for 
the I = 2-correlators must be considered as fortuitous be- 
cause the overestimation of the nonergodicity parameters 
(seen for I = 0, 1) is compensated by a generally too small 
amplitude of the I — 2-correlators. 
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FIG. 3. Diagonal critical nonergodicity parameters 
F2o,2o(?, to) as calculated from the MD simulation (symbols) 
compared with the theoretical predictions of MMCT in diag- 
onal-dipole approximation (short dashed line), in dipole ap- 
proximation (solid line) and the results obtained using the 
"full" set of MMCT-equations (long dashed line). 

The reason for the worse agreement of the I = 2 
correlators is their higher sensibility to the cutoff at 
l co = 2 which can be understood on a mathemati- 
cal level by a closer examination of the vertices which 
have been given in ref. (!(]]. Let us pick out one spe- 
cial example to illustrate this point. The vertex factor 

V 00,20,20 



(q, qi,q2] 0, 0, 0) is responsible for the coupling of 



two correlators involving I 
is of the form 

u oo!20,2oO?'9i,92;0,0,0) : 



and 1 = 2, respectively. It 



(9) 



where Uqq 1u 2 o(q, Ql, <?2! 0, 0, 0) contains Clebsch-Gordan 
coefficients of the form C(l, 2, 0; m', — to', 0) which en- 
forces 1 — 2. Therefore the contribution to the mem- 
ory function matrix caused by the vertex factor in Eq. 
@ is exact even for l co = 2. The vertex factor 
(Qi Qii 12] 0, 0, 0) instead, which describes the cou- 



^20,20,20 



pling of two correlators involving I 
tively, has the form 



2 and I 



respec- 



20,20,20 '92 ; 0, 0, 0) = (10) 

X^ M 2o M ;«,2o('?' '7i J 92;0,0,0)c/„ !2 o(gi,0) + (1 <-> 2). 

l.n 

The corresponding Clebsch-Gordan coefficient 
C(Z,2,2;m', -m',0) allows I E {0,1,2,3,4}. Thus for 
I co = 2 this vertex is not entirely taken into acount, 
although it is responsible for the coupling of I = 2- 
correlations. From this example we can see that the 
I = 2-quantities are more sensitive to the cutoff than 
those for I = 0. 

The off diagonal critical nonergodicity parameters with 
n = n' = shown in figure || are reproduced well by 
the theory. The difference between dipole-approximation 
and "full" theory are relatively small. 




q (nm ) q (nm 

FIG. 4. Off diagonal critical nonergodicity parameters 
Fia,vo(. a i m ) with ( / (' as calculated from the MD simu- 
lation (symbols) compared with the theoretical predictions 
of MMCT in dipole approximation (solid line) and the re- 
sults obtained using the "full" set of MMCT-equations (long 
dashed line). 

Figures ^ and ^| show the comparison between the new 
critical nonergodicity parameters with n ^ and/or n' ^ 
and the simulation results. Since those terms 
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FIG. 5. Diagonal critical nonergodicity parameters 
F22, 22 (a, m) with n 7^ as calculated from the MD simulation 
(symbols) compared with the results obtained using the "full" 
set of MMCT-equations (long dashed line). 
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FIG. 6. Diagonal critical nonergodicity parameters 
-p2n,2n' (? 5 m ) with ti ^ and/or n' / as calculated from 
the MD simulation (symbols) compared with the results ob- 
tained using the "full" set of MMCT-equations (long dashed 
line) . 

are neglected by the approximations they can only be 
calculated using the "full" set of MMCT-equations. The 



agreement between theory and simulation is satisfactory 
for almost all correlators. 

In the discussion of the unnormalized critical noner- 
godicity parameters we have seen that the worse perfor- 
mance of the " full" theory is partly due to the fact that 
the static input has to be taken at the "wrong" temper- 
ature. This influences the nonergodicity parameters in 
two ways: i) S(q, m) is the initial value of the correlation 
function S(g, m, t) whose limit for t — > oo are the un- 
normalized nonergodicity parameters Fi nt i> n > (q, m). ii) 
In the mode coupling approximation the static structure 
factors determine the vertices which describe the cou- 
pling between the tensorial density modes in the system. 
This influence can partly be eliminated by calculating the 
normalized critical nonergodicity parameters 



fln,Vn' (q,m) 



In.l' 



1 (<?, rn) 



n' ,l'n' 

(q,m) 



(11) 
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FIG. 7. Normalized critical nonergodicity parameters 
foo,oo(q, 0) for the centre of mass. The results from the MD 
simulation (symbols) are compared with the theoretical pre- 
dictions of MMCT in diagonal-dipole approximation (short 
dashed line), in dipole approximation (solid line) and the re- 
sults obtained using the "full" set of MMCT-equations (long 
dashed line). 

A selection of the normalized diagonal nonergodicity 
parameters is shown in Fig. (]?])- (||). Due to the normal- 
ization the variation with q is less pronounced. We ob- 
serve that the difference in amplitude between the " full" 
theory and the approximations is reduced as the initial 
value S(q,m,t = 0) of the time-dependent correlation 
functions is set to 1 at all wavevectors and all temper- 
atures. Nonetheless the dipole and diagonal-dipole ap- 
proximation still provide the better description of the 
normalized critical nonergodicity parameter. 
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FIG. 8. Diagonal normalized critical nonergodicity param- 
eters /oo,oo((Z, m) as calculated from the MD simulation (sym- 
bols) compared with the theoretical predictions of MMCT in 
diagonal-dipole approximation (short dashed line), in dipole 
approximation (solid line) and the results obtained using the 
"full" set of MMCT-equations (long dashed line). 
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FIG. 9. Diagonal normalized critical nonergodicity param- 
eters /2o,2o(<Z, m) as calculated from the MD simulation (sym- 
bols) compared with the theoretical predictions of MMCT in 
diagonal-dipole approximation (short dashed line), in dipole 
approximation (solid line) and the results obtained using the 
"full" set of MMCT-equations (long dashed line). 

In part this may be explained by the better convergence 



of the iteration for the approximation schemes, but the 
wrong critical temperature also has its influence as can be 
seen from the fact that the peak positions for the "full" 
MMCT results are still slightly shifted. 

III. SUMMARY AND CONCLUSIONS 

From our analysis we can state the following conclu- 
sions: i) The mode coupling theory for molecular liquids 
overestimates the critical temperature T c in the same 
fashion as its counterpart for simple liquids, ii) The in- 
dependence of the critical nonergodicity parameters is 
well reproduced in the vicinity of the main peaks. Sys- 
tematic differences exist for large wave vectors, in) De- 
viations between theory and simulation are partly due to 
the overestimation of the critical temperature. Thcrforc 
the structural input of the theory does not reflect prop- 
erly the structure of the liquid at the true glass tran- 
sition temperature. These deviations still exist for the 
normalized nonergodicity parameters, iv) Approxima- 
tion schemes taking into account only part of the cor- 
relators like the dipole and the diagonal-dipole approxi- 
mation can already give a reasonable description, v) The 
essential correlation functions can be selected on the basis 
of the static structure factors. I.e. for supercooled wa- 
ter the most important static correlator are those with 
n 7^ and/or n' ^ and for n = n' = those, which are 
diagonal in I and I'. 

We think that in combination with suitable approxi- 
mations for the static structure factors MMCT will be 
a useful tool in the understanding of the glass transition 
in molecular liquids. The solution of the "full" set of 
MMCT equations is still a formidable task with present 
day computers, but approximation schemes can be con- 
structed on the basis of the importance of various static 
structure factors. A reduction to the most prominent cor- 
relators and g-vectors offers the possibility to construct 
microscopically based schematic models for the descrip- 
tion of molecular liquids. 
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